Dual Nature of Improper Ferroelectricity in a Magnetoelectric Multiferroic 
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Using first principles calculations, we study the microscopic origin of ferroelectricity (FE) induced 
by magnetic order in the orthorhombic HoMnOa. We obtain the largest ferroelectric polarization 
observed in the whole class of improper magnetic ferroelectrics to date. We find that the two 
proposed mechanisms for FE in multiferroics, lattice- and electronic-based, are simultaneously active 
in this compound: a large portion of the ferroelectric polarization arises due to quantum-mechanical 
effects of electron orbital polarization, in addition to the conventional polar atomic displacements. 
An interesting mechanism for switching the magnetoelectric domains by an electric field via a 180° 
coherent rotation of Mn spins is also proposed. 

PACS numbers: 75.80.+q, 75.50.Ee, 77.80.-e, 75.47.Lx 
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Magnetoelectric materials owing their improper ferro- 
electric (FE) order to symmetry breaking magnetic struc- 
tures are drawing enormous recent interest [l], 0, H] • One 
of the fundamental problems in this area is the under- 
standing of the microscopic origin of their electric polar- 
ization. Two basic mechanisms have been proposed in 
model studies. According to one of them, magnetic or- 
dering results in the hybridization of electronic orbitals 
producing a polar charge distribution (3,0]. The other, 
more conventional approach, views the displacements of 
ions from their centrosymmetricpositions as the primary 
source of the polarization @, 0, Q. Extensive experi- 
mental studies have not been able to distinguish between 
the two possibilities due to very small values of the po- 
larization P found in this class of multiferroics, such as 
P < 0.1 ^iC/cm 2 in TbMn0 3 and TbMn 2 5 0,(1]. 

In a quest for higher P, a recent model Hamiltonian 
study concentrated on the collinear antiferromagnetic-E 
(AFM-E) spin configuration, where ferromagnetic zigzag 
spin-chains in the Mn02 planes are antiferromagneti- 
cally coupled with respect to both adjacent in-plane 
chains (see Fig. Qii) and out-of-plane stacked chains, 
as found in the orthorhombic HoMn03 and other per- 
ovskite compounds [1, 0]. The predicted polarization 
P = 0.5 — 12/xC/cm 2 was much higher than in other 
improper magnetic ferroelectrics. However, pyroelectric 
current measurements on bulk polycrystalline samples re- 
vealed ferroelectricity in HoMnC>3 with P of only less 
than 2 nC/cm 2 [l(J. 

First-principles calculations appear to be well-suited 
to address the two issues discussed above: (1) they 
can reveal the dominant mechanism of improper FE 
polarization in magnetically ordered compounds, since 
the electronic structure and lattice distortions can be 
treated simultaneously, and (2) they can solve the theory- 
experiment discrepancy for HoMnOain the AFM-E phase 
by clarifying the true value of the electric polarization. 




FIG. 1: (Color online) a) The in-plane arrangement of Mn 
and O atoms. Arrows denote the direction of spins and 
antiferromagnetically-coupled zig-zag spin-chains are high- 
lighted by pink and blue shaded areas. Structural parameters 
reported in Table [T] are shown, b) Arrows show the directions 
of the in-plane ionic displacements for Mn (left) and O (right) 
in the AFM-E phase. The thick arrows at the bottom show 
the direction of the resulting displacements the Mn and O 
sublattices and P. 



We perform simulations based on the generalized gra- 
dient approximation (GGA) [ll[ to density functional 
(DF) theory using the Vienna Ab-initio Simulation Pack- 
age (VASP) 11211 and the projector-augmented- wave pseu- 
dopotentials [l3j . Ho 4/ electrons were assumed as frozen 
in the core. The plane wave energy cut-off was set to 
500 (400) eV for the collinear (non-collinear) calcula- 
tions. The Brillouin zone sampling was performed using 
the 3x4x6 shell Q. The GGA+U calculations are 
performed by applying a Hubbard-like potential for Mn 
d states, with U ranging from to 8 eV; accordingly, we 
chose J = 0.15 U. The Berry phase approach 1(| 17 1 was 
used to calculate P, integrating over six k-point strings 
parallel to the c axis, each string containing 6 k-points. 
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Non-collincar calculations were performed according to 
Ref. [IH. Spin-orbit coupling (SOC) was neglected. As 
for the structural details, we chose the experimental lat- 
tice constants [19] for HoMnC>3 for the orthorhombic 
unit-cell (space group Pnma, a = 5.835 A, b = 7.361 A 
and c = 5.257 A) and performed atomic relaxations until 
the Hellman-Feynman forces were below 0.015 eV/A. 

Due to the small size of the Ho ionic radius, HoMn03 
shows [l9| a highly distorted perovskite structure with 
the Mn-O-Mn ac-plane angle a ~ 144°. The AFM E- 
type spin-alignment is stabilized due to this strong dis- 
tortion of the perovskite structure, as was s ugg ested by 
previous model 2fJ, 21 1, and first-principles 22 1 reports. 
With the optimized atomic positions obtained after im- 
posing the spin order, we calculated the FE polariza- 
tion using the point charge model (PCM) with nominal 
charges. The calculated Ppcm is shown in Fig. [2p as 
a function of the Hubbard parameter U. Ppcm{U) is a 
decreasing function, which is consistent with the general 
ideas of the model calculations [8]. In that model, the 
polarization appears due to the difference (a p — a ap ) be- 
tween the Mn-O-Mn angles corresponding to the bonds 
with parallel (a p ) and antiparallel (a ap ) spins, which is a 
consequence of the Hund's coupling and virtual electronic 
hopping. In turn, U characterizes the energy penalty 
paid for adding an additional electron on a Mn site. 
Thus, increasing U makes the virtual electron hopping 
less favorable, which reduces a p — a ap and, ultimately, 
P. Nevertheless, as shown in Fig. [21 the effect of even a 
very large U = 8 eV is to decrease P by less than half an 
order of magnitude. As the value of U for HoMnC>3 is not 
known from experiments, we will resort to a parameter- 
free DF treatment with U = case from now onwards, 
noting that the calculated quantities should be trusted 
at least with respect to their orders of magnitude. 

A close look at the geometrically optimized structure 
in the AFM-E phase reveals a complicated picture of dis- 
placements shown in Fig. Q}). The Mn atoms displace 
by 0.04 A with respect to the initial centrosymmetric 
Pnma structure. The a and b components of these dis- 
placements of different Mn atoms compensate each other, 
whereas the c components add up to a net displacement 
of 0.01 A per Mn atom along the negative direction of the 
c-axis. Similarly, the in-plane O atoms are displaced on 
average by 0.02 A with the resulting c-axis displacement 
of 0.01 A per atom in the positive direction. Taking into 
account relatively smaller c-axis displacements of inter- 
plane O and Ho atoms we obtain Ppcm = 2.1/^C/cm 2 . 
Therefore, our calculations with no fitting parameters in- 
dependently confirm the large electric polarization ob- 
tained in the previous model [8[. In this regard, we hope 
that the experimentally attained P could be improved 
substantially as it was accomplished for another promis- 

H0. 



ing multiferroic BiFeC>3 

To have a better understanding of the structural dis- 
tortions caused by the magnetic order, we also optimize 




FIG. 2: (Color online) a) Mn-O-Mn angles (in degrees) in the 
parallel (a p ) and antiparallel (a ap ) Mn spin configurations vs. 
the Hubbard parameter U. b) Polarization calculated by the 
point charge model and a p — a ap in the AFM-E phase. 



TABLE I: Relevant structural parameters for the AFM-A, G 
and E spin configurations: Mn-O-Mn angle (in degrees) for 
parallel (a p ) and anti-parallel (a ap ) Mn spins, large (d { and 
df) and small {d\ and d 2 s ) Mn-O bond-lengths (in A). 





a p 




d\ 


df 


dl 


dl 


AFM-A 


143.8 




2.20 


2.20 


1.93 


1.93 


AFM-G 




142.8 


2.24 


2.24 


1.90 


1.90 


AFM-E 


145.3 


141.9 


2.25 


2.18 


1.92 


1.92 



the structural parameters for the paraelectric AFM-A (all 
spins in the ac-plane are parallel) and AFM-G (all neigh- 
boring spins are antiparallel) phases and compare them 
to the AFM-E structure (see Tabic |T|. In agreement with 
the above discussion, the calculated a p for AFM-A is 
larger then a ap for AFM-G. However, when both types 
of bonds (parallel and antiparallel) are present in the 
AFM-E phase, this difference is even more pronounced, 
which explains the relatively high P. In addition, the 
[/-dependence of the angles is shown in Fig. 

The AFM-E phase shows two different kinds of AFM 
domains Q, Ei and E 2 (see the left and right insets of 
Fig. [3^, respectively) , which are expected to show oppo- 
site polarization, —P c and P c . In our calculational unit 
cell, Ei and E2 differ in the orientation of half of the spins 
(see grey highlighted regions in the central inset). Here, 
we consider a FE-AFM switching path from — P c (E\) to 
Pc (E2) via a progressive rotation of the central spins. 
According to the basic displacement-like mechanism for 
polarization, we expect P to switch from negative (in E± ) 
to positive (in E2) and to vanish when the relative orien- 
tation of the central spins with respect to the fixed spins 
is close to 90°. The 90° spin-configuration denoted as _L 
is shown in the central inset of Fig. [3ju More precisely, 
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_L is an example of a spiral magnetic structure similar to 
that in TbMn03, but commensurate with the modula- 
tion vector k = (1/2, 0, 0), and which should be FE with 
along the c-axis 0, @, 3, 26l. Based on the macro- 
scopic symmetry considerations [27j, the components of 
the polarization vector P can be expressed as follows, 

P c = Xz{c xz sin (j>-c Q cos cj>), P a = c' xz Xx sin c6, P b = 

(1) 

where \x and \z are the components of the dielectric 
susceptibility along the a and c-axes, respectively, and <fi 
is the rotation angle of the central spins. The coefficient 
Co stems from nonrelativistic interactions, while c xz and 
c' xz originate from the coupling of P to the product of 
the a and c components of the Mn spins, which has a 
relativistic origin. 

Equations (fTJ) lead to several important conclusions. 
First, for the commensurate spiral state _L {4> = n/2), 
a longitudinal component P a of the uniform polarization 
is present in addition to P c . Second, P± is finite due to 
purely relativistic effects, in agreement with previous mi- 
croscopic models [1, @]. Third, since the relatively small 
relativistic effects such as SOC are neglected in our com- 
putations, we observe that only the c component of the 
calculated P is finite for all <j> and P oc cos <fi, in excellent 
agreement with the numerical results in Fig. [3)3, which 
are discussed below. Also, in our computations P± = 0, 
and _L is taken as the reference paraelectric structure with 
centrosymmetric positions (csp). 

If FE switching is to occur, as the spin rotation pro- 
ceeds starting from E\ , the total energy is expected to in- 
crease up to a maximum corresponding to the "paraelec- 
tric" state (_L), and then to decrease again until the E 2 
minimum is reached. Indeed, this happens when we per- 
form the non-collinear calculations by varying <f> between 
and 180° with the full optimizations of the internal 
atomic coordinates for each spin configuration. The cal- 
culated total energy of the system as a function of <j> (see 
Figj3] a) clearly shows a double-well structure, with the 
depth of the well of about 8 meV/formula unit (f.u.). Al- 
though the exact magnitude of the depth of the well can 
be affected by computational details and approximations, 
we expect the feasibility of the magnetoelectric switching 
by the application of realistic electric fields. The calcu- 
lated energy barrier is, in fact, smaller than in proper FE 
BaTi0 3 (18 meV/f.u.) and PbTi0 3 (200 meV/f.u.), and 
multiferroic BaMF 4 (> 20 meV/f.u.) [Si, HI- 

The evaluation of P deserves a careful discussion, since 
it leads to an intriguing outcome. In Fig|3jD, we re- 
port the polarization evaluated by the PCM and Berry- 
phase (BP) approaches within the density-functional- 
theory (DFT-BP) along the previously mentioned switch- 
ing path. There is a marked disagreement between the 
PCM and DFT-BP approaches, therefore suggesting that 
purely quantum electronic effects are at play in determin- 
ing the final P, similar to the conventional FEs [3(J. To 
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FIG. 3: (Color online) a) Total energy as a function of the 
orientation of the central spins (see gray highlighted regions 
in the central _L spin-configuration) with respect to the spin 
of the Mn in the origin, b) Polarization calculated via the 
PCM (empty blue diamonds) and quantum-mechanically via 
DFT-BP (filled red diamonds) . The pink triangle denotes the 
values obtained via the DFT-BP approach for centrosymmet- 
ric atomic configurations as explained in the text. The lines 
are fits to P oc — cos d> with constant coefficients. 



investigate the purely electronic effects, we calculateP 
considering the atomic positions of the _L structure and 
artificially switching the spin-configuration, without re- 
laxing the lattice degrees of freedom. In this case, due to 
structural centrosymmetry, there is no contribution from 
atomic displacements. However, the calculated BP po- 
larization is found to be up to 3.5 /LtC/cm 2 (see triangles 
in Fig. IHt>)- This large contribution arises solely from the 
electronic contribution due to symmetry breaking by the 
AFM-E ordering. 

The quantum effects can be quantified further by con- 
sidering the deviations of the effective Born charges from 
their nominal valencies. In the _L case, when each O is 
coordinated to two Mn atoms with perpendicular spins, 
we obtain Z^(Mn) = 3.9 e and Z* ± (0) = -3.1 e. Along 
the AFM-FE switching path, the O ions O p and O ap 
become increasingly different and acquire different Born 
charges due to the orientation of the Mn spins to which 
they are bonded. In the extreme points corresponding 
to AFM-E1 and AFM-E2, we obtain: Z*(Mn) = 3.8 e, 
Z*(O a ) = -2.6 e and Z*{O ap ) = -3.5 e. Note, that in 
both _L and E-type spin arrangements, the Born charges 
are not extremely different from their nominal valence, 
consistent with a rather ionic nature of the chemical 
bonds. We note that the two different kinds of O ions are 
responsible for ferroelectricity in the case of centrosym- 
metric positions, where the displacement mechanism is 
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FIG. 4: (Color online) The ac-plane charge density isosurface 
plot in the energy range [-0.8:0] eV (considering the top of 
the valence band as zero of the energy scale) for fully relaxed 
positions in the AFM-E1 spin configuration. Blue (green) 
balls denote Mn atoms with up (down) spins. 



switched off. When atomic relaxation is allowed, all the 
atoms (including displaced Mn) contribute to the final 
FE polarization. The inequivalence of O a and O ap is fur- 
ther confirmed by the charge density plot for the fully- 
relaxed structure, corresponding to the energy range of 
hybridized Mn e g and O p orbitals located just below 
the valence band maximum, as shown in Fig. 3J In ad- 
dition to the expected G-like orbital ordering 0, l20j ]. 
Fig. Q] clearly demonstrates that there is a strong asym- 
metry in the charge distribution between the two O ions. 
Moreover, focusing on O ap , the charge seems to favor 
the short Mn-0 bond compared to the long bond. This 
suggests that the polar charge distribution and related 
wave-functions are due to a delicate combination of the 
Jahn- Teller effect and symmetry-breaking magnetic or- 
dering. 

In summary, our first-principles results show that, in 
the AFM-E-type HoMn03, the symmetry-induced in- 
equivalence of the in-plane Mn-O-Mn configurations for 
parallel and antiparallel spins is an efficient mechanism 
in driving a considerable ferroelectric polarization. The 
calculated total polarization of the AFM-E phase, P w 
6/iC/cm 2 is consistent with the previously obtained the- 
oretical estimates Q. In addition to the displacement 
mechanism, we find a larger but comparable contribu- 
tion arising from a purely electronic quantum effect of 
orbital polarization. The finite ferroelectric polariza- 
tion, even with a centrosymmetric atomic arrangement, 
is an unambiguous indication of a magnetism-induced 
electronic mechanism at play. Also, a magnetoelectric 
domain switching path is proposed, in which the re- 
versal of polarity of the applied electric field induces a 



180°-nip of selected spins. Although we have focused on 
the case of HoMnC>3 as example, we believe our results 
concerning the dual nature of ferroelectricity as arising 
from a symmetry breaking induced by the magnetic order 
should have a wider validity for improper magnetic ferro- 
electrics. Our findings suggest that the interpretation of 
experiments, as well as model calculations, should take 
into account both the lattice and electronic mechanisms 
of improper ferroelectricity in multiferroics. 
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Auxiliary Material 



Here we give a brief derivation of selected results of the Landau theory theory of phase transitions [3l( applied to 
the AFM-E order and its coupling to ferroelectricity. Here we work in the Pnma (#62) space group setting, which is 
related to Pbnm used in Ref. [S| via a different choice of the orthorhombic axes[32|. We also place the origin in the Mn 
position. The order parameter of the magnetic ordering (Ei,E2) is written as a linear combination of Mn spins S^, 
where i — 1, . . . , 8 denotes the Mn atoms occupying the following positions in the magnetic unit cell, corresponding 
to the modulation vector k = (1/2, 0, 0): 



1:(000), 2:(-0-), 3 : (100), 
The order parameter is obtained asQ 



,3 1 



1 1 1. 



'2 2 



4:^0-), 5: (0-0), 6 : (---), 

S7 — Sg, 

s 7 + s 8 . 



7:(y)), 



,3 1 1 
'2 22 



Ei — Si — S2 — S3 + S4 — S5 + Sg 
E2 = Si + S2 — S3 — S4 — Sg — Sg 



(2) 



The components of Ei and E2 transform according to the rules given in Table HT1 where the transformation rules for 
electric polarization P are also summarized. Here, we analyze the possible magnetoelectric coupling terms of the form 
P ■ E 2 in the thermodynamic potential: 



F me = c P z (El--E 2 2 ) 



■ yy 



y2 



E. 



1/2 \ 



Z P Z {E? - E?) + c xz P z {EfE* - E x 2 El 



+d xy P y {E v 1 E% ~ E\Ef) + d zy P y (E\El - E V 2 E Z 2 ) 



c' xz P x (E*E% - EfE?) 

(3) 



Clearly, only the first term has a nonrelativistic origin since (Ei — E 2 ) is a linear combination of the terms (Si • S&). 

Considering Eq. ((3]) together with the usual dielectric energy Fd = l/2(x x ~ 1 P 2 + Xy X Py + Xz 1 Pz)> we obtain that 
P y = if E\ — E 2 = 0, which always holds in our case. The other two components of P are given as 



P x = -X x c' xz {E* 2 E% - EfEf), P z = -x.ME? - Ei) 



(£f - E?) + c xz (EfEI - E$Ef) 



(4) 



As is proposed in the main text, the switching of the magnetoelectric domain is accomplished by rotating the spins 
S2, S4, Sg, and Ss in the xz-plane. Thus, we set 



Si 



Using Eqs. ([2]) and (J4J) , it is easy to verify that 
respectively. Furthermore, we obtain, 



5(1,0,0), -S 2 
= and 6 - 



- S 4 = Sg = — S 8 = 5(cos0, 0,sin0), 
7r correspond to the domains Ef(—P z ) and E 2 ( 



-P, 



P x = Xx^c' xz sin4>, P z = x»(32S <, c x ,Bin0 -64^00 cob «^), 

which are equivalent to the expressions given in the main text with the numerical constants and S 2 adsorbed in the 
coefficients. 

As a by-product of the present analysis we also find that the following term is allowed in the thermodynamic 
potential: 



F L = A J(E*d x E% - E$d x E?)dV. 
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TABLE II: Matrices of the generators of space group Pnma in the representations spanned by Ei, E2, and P. The space group 
elements are denoted (r\hkl), where r is the identity operation 1, two-fold rotation 2 a ,c, inversion I, or time reversal 1' followed 
by the translation r = ha. + kh + Ic. 





(2.|i0|) 

V " 1 2 A t 


(2 y |0i0) 


(I|000) 


(1|010) 


(l'|000) 


Ef 


1 


1 


1 


-1 


-1 


E% 


-1 


1 


1 


-1 


-1 


El 


1 


-1 


1 


-1 


-1 


El 


-1 


-1 


1 


-1 


-1 


El 


-1 


1 


1 


-1 


-1 


El 


1 


1 


1 


-1 


-1 


Px 


-1 


-1 


-1 


1 


1 


Py 


-1 


1 


-1 


1 


1 




1 


-1 


-1 


1 


1 



This Lifshitz [3lL l33l ] term leads to the following typical series of phase transition when the temperature is low- 
ered: Para-phase — > Incommensurate phase — > "Locked- in" commensurate phase [33j|. which was indeed observed 
in HoMn0 3 @. It can be shown[33( that both Ef and #f are modulated with k(T) = (1/2 + 5k(T), 0,0) in the 
AFM-Incommensurate phase. Taking into account the magnetoelectric interactions ([3]), this leads to a vanishing 
macroscopic polarization in the incommensurate phase in a greement with the experiment jlO|. 



